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Abstract. The Casimir effect, reflecting quantum vacuum fluctuations in the 
electromagnetic field in a region with material boundaries, has been studied both 
theoretically and experimentally since 1948. The forces between dielectric and metallic 
surfaces both plane and curved have been measured at the 10 to 1 percent level 
in a variety of room-temperature experiments, and remarkable agreement with the 
zero-temperature theory has been achieved. In fitting the data various corrections 
due to surface roughness, patch potentials, curvature, and temperature have been 
incorporated. It is the latter that is the subject of the present article. We point 
out that, in fact, no temperature dependence has yet been detected, and that the 
experimental situation is still too fluid to permit conclusions about thermal corrections 
to the Casimir effect. Theoretically, there are subtle issues concerning thermodynamics 
and electrodynamics which have resulted in disparate predictions concerning the nature 
of these corrections. However, a general consensus has seemed to emerge that suggests 
that the temperature correction to the Casimir effect is relatively large, and should be 
observable in future experiments involving surfaces separated at the few micrometer 
scale. 
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1. Introduction 

About the same time that Schwinger and Feynman were inventing renormalized quan- 
tum electrodynamics, Casimir discovered that quantum electrodynamic fluctuations re- 
sulted in macroscopic forces between conductors and dielectrics pQ. The theory was a 
natural outgrowth of the Casimir-Polder theory of the retarded dispersion force between 
molecules [2]. The general theory for the forces between parallel dielectrics was worked 
out by Lifshitz and collaborators [3], who also included temperature corrections, which 
were considered further by Sauer jl] and Mehra [3]. Some years later, the whole theory 
was rederived by Schwinger, DeRaad, and Milton [B]. 

The early experiments on Casimir forces were rather inconclusive - for a review 
see [7]. However, the corresponding Lifshitz theory was verified rather impressively 
by Sabisky and Anderson [Hj, so there could hardly be any doubt of the validity of 
the essential ideas. Starting about a decade ago, modern experiments by Lamoreaux 
CHI Ell E!, Mohideen and collaborators [TJ Q3J US!, and by Erdeth PU brought 
the experimental measurement of the Casimir force between curved metal surfaces 
(mapped to the plane geometry by the proximity approximation ^3QB]) into the percent 
accuracy region. (Exact results have now apparently rendered the use of the proximity 
approximation, which cannot be extended beyond leading order, unnecessary. See, for 
example, ^01 1201 E]-) Application of such Casimir forces to nanoelectromechanical 
devices have been suggested by experiments at Bell Labs and Harvard |221 I2H1 E3] Only 
one experiment so far, of limited accuracy (~ 15%), has employed parallel plates [2.5] • 
The difficulty of maintaining parallelism in that geometry limits the accuracy of the 
experiment, but the forces are much larger than those between a sphere and a plate, 
so the forces can, in principle, be determined at much larger separations. Proposals to 
perform measurements of the force between a cylinder and a plane [2H] and between 
eccentric cylinders j2Ij have advantages because the forces are stronger than between 
a sphere and a plane, yet the difficulties in assuring parallelism are not so severe as 
with two plane surfaces. The most precise experiments so far, based on both static and 
dynamical procedures between a plate and a spherical surface, have been performed at 
Purdue j2HH23E01, where the accuracy is claimed to be better than 1% at separations 
down to less than 100 nm. 

All present experiments agree well with the zero-temperature Casimir theory when 
surface roughness and finite conductivity corrections are included [HU E2j- The issue 
about which controversy has recently erupted is the temperature dependence. (For 
recent statements of both sides of the controversy, see jHOl E23 EH EEl EE]-) All 
experiments reported to date have been conducted at room temperature, so there is 
no direct evidence for or against any particular model of the temperature dependence. 
Indirect evidence for this dependence has been inferred based on the nonzero shift in 
the theoretical Casimir force between the surfaces due to the difference between the 
force at zero temperature and at 300 K. Surprisingly, this temperature shift is not so 
straightforwardly computed as one would have at first suspected. 
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It is the purpose of the present paper to frame the question of the temperature 
dependence of the Casimir force in the context of the history of the subject and 
the present experimental constraints, as well as to point out ways of reconciling the 
ambiguities both from the theoretical and experimental sides. In the following section 
we review the standard approach given in [Hj for both dielectric and metal surfaces. 
Then, in section |3] we give the arguments why the transverse electric (TE) zero mode 
should not be included, and how this impacts the temperature dependence of the 
force, and the resulting impact on the free energy and entropy. Other theoretical 
arguments for and against this point of view are discussed in section HJ The status 
of the experimental situation, and the possibility of dedicated experiments to search for 
the temperature dependence of the Casimir effect, will be reviewed in section |S] Some 
new calculations are presented in section |H1 in the hope of providing signatures to help 
resolve the controversy. Finally, concluding remarks are offered in section |3 



2. Conventional temperature approach 

The zero-temperature Casimir effect between parallel conducting plates, or between 
parallel dielectrics, is very well understood, and is not controversial. The formula for 
the latter, which includes the former as a singular limit, may be derived by a multitude 
of formalisms, which will not be reviewed here [HJ EHJ EUl EEl EHI • For a system of parallel 
dielectric media, characterized by a permittivity 



e(z) — I £ 3 , < z < a, (2.1) 




where the various permittivities are functions of frequency, the Lifshitz force per unit 
area on one of the surfaces is at zero temperature 

1 poo poo 

pT=0 = -^J o dC ^ dkln^d-' + d'- 1 ), (2.2) 
where ( is the imaginary frequency, ( = —iu, and the longitudinal wavenumber is 



/;„ Jkl + C 2 ei(K), (2.3) 

while the transverse electric (TE) and transverse magnetic (TM) Green's functions are 
characterized by the denominators 

d= V K3a -l, d' = , ] , V K3a -l, (2.4) 

K3 — K\ K3 — K2 K 3 — K\ K 3 — K 2 

respectively, where k[ = K^jsi. 

The attractive Casimir pressure between parallel perfectly conducting planes 
separated by a vacuum space of thickness a is obtained by setting e 1>2 — > 00 and e 3 = 1. 
In that case the TE and TM contributions are equal, and we have 

Pc = ^ L dC L ^e^I = ~-2 / d C^T37 = -240^' (2 ' 5) 
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which is Casimir's celebrated result PQ. 

The controversy surrounds the question of how to incorporate thermal corrections 
into the latter result. At first glance, the procedure to do this seems straightforward. It 
is well-known that thermal Green's functions must be periodic in imaginary time, with 
period (3 = 1/T [10] • This implies a Fourier series decomposition, rather than a Fourier 
transform, where in place of the imaginary frequency integral in ()2.5|) we have a sum 
over Matsubara frequencies 

C = (2-6) 
that is, the replacement 

M m=0 

the prime being an instruction to count the m = term in the sum with half weight. 
This prescription leads to the following formula for the Casimir pressure between perfect 
conductors at temperature T, 

"OO 



-1 v P OO 1 



where 



r m=0 

From this it is straightforward to find the high and low temperature limits, 

pT ~-i^ c(3) -2^( i+i 4) e "' " <<4 ™> (2io » ) 



p 



T 



71 



2 



16 a* _ 240 a p/a 
+ 3 ^ tt f3 



(3 > 4tto. (2.106) 



240a 4 

These are the results found by Sauer |lj and Mehra and by Lifshitz [H]. The two 
limits are connected by the duality symmetry found by Brown and Maclay The 
pressure may be obtained by differentiating the free energy, 

p = ~l F ' < 2 - n > 

which takes the following form for low temperature (now omiting the exponentially small 
terms) 

^ 2 ((3)^3 , ^ 



T i + t\ aT < 1, (2.12) 



2 



720a 3 2tt 45 
from which the entropy follows, 

s ~-W F ~ 3j ¥ T2 -lK T3a < aT<<1 - (2 ' 13) 

which vanishes as T goes to zero, in accordance with the third law of thermodynamics, 
the Nernst heat theorem. 
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3. Exclusion of TE zero mode 

However, there is something peculiar about the procedure adopted above for a perfect 
metal. (This seems first to have been appreciated by Bostrom and Sernelius [42J.) It has 
to do with the transverse electric mode of zero frequency, which we shall refer to as the 
TE zero mode. If we examine the zero frequency behavior of the reflection coefficients 
for a dielectric appearing in (|2.4j) . we see that providing C^eiiQ — > as ( — > 0, the 
longitudinal wavenumber C — > 0, and hence d — > oo as C — > 0. This means 

that there is no TE zero mode for a dielectric. This statement does not appear to be 
controversial [13]. However, if a metal is modeled as the e — > oo limit of a dielectric, 
the same conclusion would apply. Because that would spoil the concordance with the 
third law noted in the previous section, the prescription was promulgated in [B] that 
the e — > oo limit be taken before the C — > limit. But, of course, a real metal is not 
described by such a mathematical limit, so we must examine the physics carefully. 
A simple model for the dielectric function is the plasma dispersion relation, 

e(u>) = \-% (3.1) 

where u p is the plasma frequency. For this dispersion relation, the condition ( 2 e(i() — > 
fails to hold as £ — > 0, and the idealized prescription result, namely the contribution of 
the TE zero mode, holds. However, real metals are not well described by this dispersion 
relation. Rather, the Drude model, 

E{K) = 1+ (iky (3 - 2) 

where the relaxation frequency 7 represents dissipation, very accurately fits optical 
experimental data for the permittivity for ( < 2 x 10 15 rad/s ^UIS]- F° r example, for 
gold, appropriate values of the parameters are [IB] 

to p = 9.03 eV, 7 = 0.0345 eV. (3.3) 

In this case, the arguments given above for the exclusion of the TE zero mode apply. 

The arguments are a bit subtle EH] , so we review and extend them here. Let 
us write the Lifshitz formula at finite temperature in the form 

x poo 00 rj 

P T = J2'f™= dm/M-^T^/^O), (3.4) 

m=0 ^° k=0 ^ '' 

where the second equality uses the Euler-Maclaurin sum formula, in terms of 

1 f°° 

f(m) = J Q dk 2 ± K (Cm) K 1 + d'- 1 ) , (3.5) 

according to (|2.2|) and (|2.7|) . where we assume that vacuum separates the two plates 
so K-s(Cm) = «(Cm) — v^I + Cm- Here the denominators ()2.4|) are functions of Cm- By 
changing the integration variable from m to Cm we immediately see that the integral 
term in the Euler-Maclaurin sum formula corresponds precisely to the zero-temperature 
result El). 
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One must, however, be careful in computing the low temperature corrections to 
this. One cannot directly expand the denominator d in powers of ( because the k± 
integral in (j3.5j) ranges down to zero. Let us rewrite the TE term there as follows: 



./ 



(TE). 



m = — 



1 

7T/3 



dn K 



2rmr//3 



1 + Vl + C(e(iC m )-l)/« 2 ' 
1- Vl + C(e(iCm)-l)/« 2 



-i 



1 



Evidently, for the Drude model, or more generally, whenever 
BmC a [e(iC)-l] = 0, 



(3.6) 



(3.7) 



/( TE )(0) = 0. However, it is important to appreciate the physical discontinuity between 
m = and m = 1 for room temperature. At 300 K, while Co — 0, Ci = 27rT = 0.16 eV, 
large compared the relaxation frequency 7. Therefore, for m > 0, 

2 -1 -1 



/ 



(TE). 



m) ~ — 



1 



dn K 



1 + u 2 Jk 2 + 1 



2rea 



1 

'tt/3 



dn K 



^2na 



- 1 



(3. 



provided the significant values of ( m and k are small compared the the plasma frequency 
Up. This is just the ideal metal result contained in (|2.8jl . Insofar as this is accurate, 
this expression yields the low- and high-temperature corrections seen in (j2.10fcj) . (j2.10aj) . 
However, there is now a discontinuity in the function / < - TE - ) . As ( rn — > 0, 



f {TE \m) - 



1 



dn- 



C(3) 



,2/ta 



(3.9) 



e™* - 1 4?r/3a 3 : 

rather than zero. This implies an additional linear term in the pressure at low 
temperatures: 

C(3) 



P 



T=0 



+ 



aT < 1. 



ina- 



(3.10) 



Exclusion of the TE zero mode will also reduce the linear temperature dependence 
expected at high temperatures, 

C(3) 



P 1 



aT > 1, 



(3.11) 



one-half the usual ideal metal result seen in (|2.10ojl - and this is indeed predicted in 
numerical results [see figure 4 of |38j for a > 5 /im, for example.] 

Most experiments are carried out between a sphere (of radius R) and a plane. In 
this circumstance, if R ^> a, a being the separation between the sphere and the plate 
at the closest point, the force may be obtained from the proximity force approximation, 

J r = 2nRF(a), (3.12) 
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F(a) being the free energy for the case of parallel plates separated by a distance a. Thus 
in the idealized description, the low temperature dependence including our linear term 
is 



-^((3)aT + _C(3)(aT) 3 - 16(oT)< 



360a 3 - ■ , aT«l. (3.13) 

Since this conversion is trivial, in the following we will restrict attention to the 
straightforward parallel plate situation. 

These results are only approximate, because they assume the metal is ideal except 
for the exclusion of the TE zero. Elsewhere, we have referred to this model as the 
Modified Ideal Metal (MIM) model jSHHSHl- Evidently, for sufficiently low temperatures 
the approximation used here, that C\ ^> 7, breaks down, the function f(m) becomes 
continuous, and the linear term disappears. Indeed, numerical calculations based on real 
optical data for the permittivity show this transition. An example of such a calculation 
is presented in section |B] There, in figure [TJ we see a negative slope in the quantity 
P/ Pc as a function of the plate separation a in the region between 1 and 2 micrometers. 
This slope is approximately — 0.1//xm. Here P is the pressure between the plates at 300 
K, while Pc is the ideal Casimir pressure (|2.5j) . If we compare this to our approximate 
prediction ()3.10jh 

p t=30ok 30 C(3) a a . . 

~p — « 1 - tt^^t— = 1 - °- 15 — > 3 - 14 

Pc i.bl 7r /xm /im 

the slope and intercepts agree at the 20% level. Accurate numerical results between real 
metal plates and sphere are given in [3*3] . 

Because this linear behavior does not persist at arbitrarily small temperatures, it 
is clear that the conflict with the third law anticipated in the arguments in the previous 
section do not apply. In fact, as we shall now see, the entropy does go to zero at zero 
temperature. 

4. Arguments in favor and against the TE zero mode 

As noted above, there are strong thermodynamic and electrodynamic arguments in favor 
of the exclusion of the TE zero mode. Essentially, the point is that a realistic physical 
system can have only one state of lowest energy. Electrodynamically, one can start from 
the Kramers- Kronig relation that relates the real and imaginary part of the permittivity, 
required by causality, which can be written in the form of a dispersion relation for the 
electric susceptibility |Hj 

xM= £r^_£M (4 „ 

An J lu' 2 — (u + \e) 2 

If the spectral function p(u') > is nonsingular at the origin, it is easily seen that 
cij 2 x(u;) — > as ui — > 0, which as shown in the previous section implies the absence of 
the TE zero mode. Conversely, p(u>') must have a 5-function singularity at the origin 
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to negate this conclusion. This would seem implausible for any but an overly idealized 
model. In contrast, in the Drude model 

It has been objected that rather than employing bulk permittivities as done in the 
usual expression for the Lifshitz formula, one should use surface impedances instead 
[IHJ HHJ EU EH] . Indeed this may be done, but it leads to identical results. The surface 
impedance merely expresses the linear relation between tangential components of the 
electric and magnetic fields at the interface between the two media, 

E ± = Z(u, kjjBi. x n, n x E ± = Z(u, k±)B ± , (4.3) 

where n is the normal to the interface at the point in question. From Maxwell's equations 
we deduce |17J EH] for the reflection coefficient for the TE modes 

r TE = -^|^, K 2 = ( 2 + kl, (4.4) 
and the surface impedance is§ 

c 



v/C 2 KiC)-i] + K 2 ' 



(4.5) 



From this reflection coefficient the Lifshitz formula is constructed according to d = 
(r TE ) e 2Ka — 1. Evidently the resultant expression for the Lifshitz pressure coincides 
with that found from the permittivity, seen for example in ()3.6|) . This coincidence 
has been well recognized by previous authors [HU E2| • The reason why Mostepanenko 
and collaborators obtain a different result is that they omit the transverse momentum 
dependence in (J4.5)) and thereby argue that at zero frequency Z vanishes, 

Z^--^~A/C, (4-6) 
which is the content of the normal skin effect formula 



where a is the conductivity. [These two formulas are seen to be identical if we replace 
oj = i( and recognize that 7 = uz,/ (Anta).] These formulas apply when we have the 
restriction appropriate to real photons k\< uj 2 . However, no such mass-shell condition 
applies to the virtual or evanescent photons involved in the thermal Casimir effect. The 
same sort of error seems to be made by Torgerson and Lamoreaux [SHI 02], and by 
Bimonte EH] • 

As noted above, use of the plasma model in the reflection coefficients would lead 
to the conventional temperature dependence, but this dispersion relation is inconsistent 
with real data. However, it has been argued that in the ideal Bloch-Griineisen model 

§ Here we have assumed that the permittivity is independent of transverse momentum. In principle 
this is incorrect, although optical data suggest that the transverse momentum dependence of e is rather 
small. See also |50| . 
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[57] the relaxation parameter goes to zero at zero temperature. However, real metals 
exhibit scattering by impurities; in any case, at sufficiently low temperatures the residual 
value of the relaxation parameter does not play a role, as the frequency characteristic 
of the anomalous skin effect becomes dominant jSH] • Moreover, the authors of [3UJ EH] 
also extrapolate the plasma formula from the infrared region down to zero frequency, 
whereas in fact frequencies very small compared to the frequency corresponding to the 
separation distance play a dominant role in the temperature dependence |58j . Finally, 
we emphasize that all present experiments are carried out at room temperature, where 
the known room temperature data are relevant. 

The principal reason for the theoretical controversy has to do with the purported 
violation of the third law of thermodynamics if the TE zero mode is not included. If 
ideal metal reflection coefficients are used otherwise (the MIM model) such a violation 
indeed occurs, because the free energy per unit area for small temperature then behaves 
like 

C(3) 
167ra 2 

However, we and others have shown I39| lo'Hl loTl] that for real metals, the free energy 
per area has a vanishing slope at the origin. Indeed, in the Drude model we have 

F = F + T 2 -^(21n2-l), (4.9) 
487 

for sufficiently low temperatures. There is, however, an intermediate range of 
temperatures where it is expected that the entropy is negative. We do not believe that 
this presents a thermodynamic difficulty, and reflects the fact that the electrodynamic 
fluctuations being considered represent only part of the complete physical system 
(SB EH IHHj , although this is not a universal opinion j^E] • (See also the further remarks 
in section [7J) New calculations are underway, showing explicitly the zero slope of the 
curve for the free energy near T = 0, thus corresponding to zero entropy 61 . 

Two other recent papers also lend support to our point of view. Jancovici and 
Samaj and Buenzli and Martin [S3] have examined the Casimir force between ideal- 
conductor walls with emphasis on the high-temperature limit. Not surprisingly, ideal 
inert boundary conditions are shown to be inadequate, and fluctuations within the walls, 
modeled by the classical Debye-Hiickel theory, determine the high temperature behavior. 
The linear in temperature behavior of the Casimir force is found to be reduced by a 
factor of two from the behavior predicted by an ideal metal, just as in (j3.11|) . This is 
precisely the signal of the omission of the m = TE mode. Thus, it is very hard to see 
how the corresponding modification of the low-temperature behavior can be avoided. 

Further support for our conclusions can be found in the recent paper of Sernelius 
[HI], who calculates the van der Waals-Casimir force between gold plates using the 
Lindhard or random phase approximation dielectric function. The central theme of his 
work is to describe the thermal Casimir effect in terms of spatial dispersion. Physically, 
spatially nonlocal effects play a role at low frequencies because charge carriers can 
move freely over large distances. Deviations from standard local electromagnetic theory 
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can then be expected. Spatial dispersion implies that the standard Fresnel equations 
no longer apply. Moreover, because of lack of experimental data one has to rely on 
theoretical model-dielectric functions. 

Sernelius finds that for large separations the force is one-half that of the ideal 
metal, just as in the calculation in Refs. jHSEB]- This agreement is not quite trivial; 
it means that the thermal Casimir effect can be explained in two, apparently unrelated 
ways: One way is to include spatial dispersion in the formalism from the beginning, 
omitting dissipation. The other way is the conventional one, namely to describe the 
thermal effect in terms of dissipation alone, by introducing the relaxation frequency 
7 as we have done above. Sernelius shows that, for arbitrary separation between the 
plates, the spatial-dispersion results nearly exactly coincide with the local dissipation- 
based results jHHl H2] • 

5. Experimental constraints 

We have marshaled theoretical arguments that seem to us quite overwhelming in favor 
of the absence of the TE zero mode in the temperature dependence of the Casimir force 
between real metal plates, which seem to imply unambiguously that there should be large 
(~ 15%) thermal corrections to the Casimir force at separations of order 1 micrometer. 
New detailed calculations based on this theory, and using real optical data for aluminum, 
are discussed in the following section. The difficulty is that, experimentally, it is not 
easy to perform Casimir force measurements at other than room temperature, so current 
constraints on the theory all come from room temperature experiments. Then all one can 
do is compare the theory at room temperature with the experimental results, which must 
be corrected for a variety of effects, such as surface roughness, finite conductivity, and 
patch potentials. A deviation between the corrected zero temperature theory and the 
room temperature observations then is taken as a measure of the temperature correction. 

The temperature correction is evidently relatively largest at the largest separations, 
where, unfortunately, the total Casimir force is weakest. Lamoreaux's early experiments 
[Oj were conducted at the 1 /zm scale, so if they were accurate to 10% they would have 
seen the effect our theory predicts, but probably, in spite of Lamoreaux's assertion, 
they were not so accurate, because few essential corrections were included jHHl- The 
experiments of Mohideen et al fBl El CHj were much more accurate, but because they 
were conducted at much smaller distances, even our rather large temperature correction 
would have remained inaccessible. It is the most recent experiments of the Purdue 
group [23 EH] that claim the extraordinarily high precision to be able to see our effect 
at distances as small as 100 nm. Indeed, they see no deviation from the corrected 
zero-temperature Lifshitz theory using optical permittivities, and hence assert that our 
theory is decisively ruled out. The effect we predict for the temperature correction is 
only 1.5% at a distance of 160 nm [HH], so the measurement must be performed at the 
1% level to see the effect there. (For the usually employed sphere-plate configuration, 
AjF/jF ps 2.5% at a = 160 nm.) Although they claim this degree of accuracy, it is 
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doubtful that they have achieved it, because, for example, to achieve 1% accuracy, the 
separation would have to be determined to better than 0.3%, or 0.5 nm at a = 160 nm. 
Since the roughness in the surfaces involved is much larger than this (see also ""Hj), and 
other corrections (such as the fact that the metallic surfaces are actually thin films, and 
the effects of surface plasmons |""E1 EH]) have not been included, we have reason to be 
skeptical of such claims [70] . 

In any case, it would seem imperative to perform experiments at different 
temperatures in order to provide evidence for or against temperature dependence of 
Casimir forces. We understand such experiments are in progress. We encourage 
experimentalists to redouble their efforts to determine the presence or absence of such 
an effect in an unbiased manner, for the issues involved touch at the heart of our 
fundamental theoretical understanding of electrodynamics, statistical mechanics, and 
quantum field theory. 



6. New calculations 



To aid in the experimental disentanglement of this effect, we have carried out new 
calculations of the Casimir force between two infinite half-spaces made of aluminum, 
separated by a vacuum space of width a. (Other recent calculations appear in 
[3*51 rfT| HB*j.) The results are shown in figure ["3 (Figures [fl"""| are taken from the 
Master's Thesis of SAE jZ2j-) Formulas made use of are read off from ()2.2j) and are now 
given in usual dimensional units: 

fr poo poo / a2 -2n a a2 -2k <j \ 

P T _ ( a ) = -—\ dC / dk ± k ± n ( — + — ] , (6.1a) 

p r>0 f a ) = -2iLy' dk ± k ±Ko At a e 9 6 9 + A ™ e - ), (6.i6) 

7T ^ ./n V 1 - A 2 E e~ 2 *oa i _ A^j^e - ° a J ' 



m=0 



where a is the separation between plates (of infinite thickness), and 

Ate = —A (6.2a) 
k - e(i()n 

A T m = — : — TTTT j (6.26) 
k + e(i£)K 



K = X lki + s(i()^, (6.2c) 



2 



Here, e(iC) = e (iC)/ e o is the usual permittivity relative to the vacuum. In the case 
of finite temperatures, ( = ( m = 2nmkBT/h, and a standard Lifshitz substitution of 
integration variables was made during calculations (see for example (3.2b) of |HE])- The 
results are plotted relative to the standard Casimir pressure Pq in (|2.5|) . Calculations 
have been carried to a relative accuracy of better than 10~ 4 . Even at T = there are 
large deviations from the ideal Casimir result at all distance scales. 
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Casimir force between two Al half-spaces 
relative to ideal result 




Figure 1. Temperature dependence of the Casimir force between aluminum plates. 



To illustrate the contributions of the TE and TM modes, figures EHU depict the TE 
and TM integrands of a Casimir pressure expression of the type 

poo poo 

Pt=o = / dC / dA; ± [J TE (iC, k ± ) + J TM (iC, (6.3) 
Jo Jo 

It is clear that the TE term in the integrand falls off rapidly to zero as £ —>■ whereas 
the TM term remains finite. The relative contributions to the pressure by the TE and 
TM modes are illustrated in figures Evidently, the contribution of the TE mode 
rapidly decreases with increasing temperature and increasing plate separation. 

6.1. Results for a five-layer model 

Because many experiments have been carried out with a conducting surface between 
parallel capacitor plates it is useful to consider the five layer geometry which has been 
treated repeatedly by Tomas in recent years 73 j . The geometry is defined by figure 
All three slabs are assumed to be made of aluminum. We assume one intermediate plate 
of width b, immersed in a cavity of total width c = a + a' + b. Between the central plate 
and the two outer semi-infinite media we assume a vacuum (e g = 1). The quantity h is 
defined as h = c — b. The quantity 5 is the deviation of the center of the plate from the 
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Figure 2. TE part of the integrand in 



TM-polarised part of integrand 
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Figure 3. TM part of the integrand in 1)6. 3fl . 
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Figure 5. TE and TM contributions to the pressure shown in figure ^ 
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Figure 6. Ratio of the TE and TM mode pressures contributing to the pressure shown 
in figure ^ 
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Figure 7. The 5-zone geometry. Here we have set s g = 1. 
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Force on centre slab in 5-zone system 



Slab: Aluminium 
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Gap material: vacuum 
Cavity width: 3|Am 
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Figure 8. Pressure on the plate at a distance 6 from the cavity center. If S < one 
gets the antisymmetric prolongation about the position 6 = 0. 



midline of the cavity. The Casimir pressure at zero and finite temperature is given by 

TM 



P T=0 (6; b,c) = — dC / dk ± *JL5 *> b > c 

*K JO JO „_ TF 



JHo(*;6,c) = — V' 



q=TE 
TM 



7T 



(6.4a) 
(6.45) 



/CO J- 
dk L ^2 Iq(Km,k±;6,b,c). 
iii q=TE 

Here the integrand is 

J ff (iC, fe±; 5, 6, c) = k ± K 2A lq A 2q (l - e - 2K2b )e- K ° h sinh 2k 5 
x [-A^e _2Ka6 + 1 - A 2 lq e- 2Koh (e- 2K2b - A 2 2q ) - 2A lg A 2g (l - e - 2 ^ h ) e -^ h C osh2/€ (5]" 1 , 

(6.5) 

The Casimir force is positive for positive 5, and is antisymmetric around the cavity 
center 6 = 0. The index q in the A's in the formulas runs over the polarizations TE and 
TM, which are given by equations (|6.2qjl . (j6.26|) . These results are equivalent to those 
found earlier by Tomas |74| . but presented in a more illustrative form. 

Numerically, we choose c = 3 /im and b = 500 nm. The pressure P(S), calculated 
from 5 = to 5 = (c — b)/2 — 50 nm, is shown in figure |H1 All calculations are done 
with an accuracy of better that 10~ 4 in the final result, which should be sufficient for 
practical purposes. Figure M shows the pressure relative to Casimir's result for ideal 
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Force on centre slab in 5-zone system 
relative to ideal result 




Figure 9. Pressure on the plate at a distance 5 from the cavity center relative 
to Casimir's result (|6.6|l for ideal conductors. For S < one gets the symmetric 
prolongation of the curve about the position 5 = 0. 



conductors, 

_7T 2 kc ( 1 1 

° ~ ~ "240" \ {h/2 + 5y ~ {h/2 - 5Y 
As for the finite temperature calculations, it has been checked that all terms in the 
sum (except for the zero mode) lie within the frequency domain covered by Lambrecht's 
data. There is thus no extra assumption made in the calculation, such as the property 
e( 2 /c 2 — > as ( — > 0, following from the Drude relation, except to ensure that there is 
no contribution from the TE zero mode. In the T = case, the Drude relation is used 
for low frequencies. The contribution to the force coming from frequencies outside the 
Lambrecht region is very small, smaller than the accuracy of the calculation. Again, 
we see large deviations from the ideal Casimir result at all temperatures, as well as 
relatively large temperature corrections, which we hope will be readily detectable in 
future experiments. 

It would be very interesting to conduct experiments which are sensitive to the 5- 
zone geometry we have considered, namely a slab in a cavity. In particular, one could 
study the case of a slab oscillating about the center of a cavity, detecting how the 
oscillation frequency varies with temperature. Using a cavity of width 2-4 micrometers, 
say, one finds that the temperature corrections are largest for slab-wall separations on 
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the order of 1 micrometer. 
7. Conclusions 

We have shown how the Casimir pressure between parallel plates can be calculated 
with the inclusion, as well as with the exclusion, of the TE zero mode. We emphasized 
(sectionHJ) that there are strong thermodynamic and electrodynamic arguments in favour 
of the latter option. As is known in general, instead of employing bulk permittivities, 
it is possible to use surface impedances instead. We pointed out that in such a case 
it becomes necessary to take into account also the transverse momentum dependence 
in the expression (j4.5)) for the surface impedance. Otherwise, if one leaves out the 
transverse momentum, one will obtain erroneous results as has often been the case in 

the literature ^DIB3|EHIH3E3IE3IE51EE!. 

We advocate the use of Drude's dispersion relation ()3.2|) throughout. The 
alternative plasma relation ()3.1|) is inconsistent with real dispersive data. One peculiar 
effect arising from use of the Drude relation is the appearance of negative Casimir 
entropy in a finite frequency interval. The physical reason for this kind of behaviour is 
the following: It reflects the fact that we are dealing with only a part of the complete 
physical system. The effect may appear counterintuitive, but is not so uncommon 
in physics after all. One cannot apply usual thermodynamic restrictions such as 
positiveness of entropy to a "subsystem" formed by the induced interaction part of 
the free energy of the full system. This issue was discussed in detail in a previous paper 
( 38]). In section IV of that paper we introduced, as an illustrative mechanical model, a 
system of two harmonic oscillators interacting via a third one. Such oscillators represent 
a simplified picture of two parallel plates interacting via the electromagnetic field. We 
calculated the classical as well as the quantum free energy of this mechanical system, 
and found there to be a finite temperature interval for which the interaction free energy 
increases with increasing temperature, thus leading to a negative interaction entropy 
term S = —dF/dT. In this way a mechanical analogy with the Casimir interaction 
energy (and corresponding interaction entropy) was demonstrated. 

The new developments of our paper are the establishment of the complete 
mathematical formalism necessary to calculate the temperature correction to the 
Casimir force. We assumed infinite, parallel, plates, with no surface roughness included, 
and assumed in all calculations the tabulated dispersive data for bulk materials. Then, 
we showed calculated results in figures ^El which are all new. They all assume 
aluminum plates. We gave results for three-layer geometry, and also some for a five-layer 
geometry. Our general suggestion for experimentally measuring the influence from finite 
temperature is to keep the geometry of the setup as constant as possible (by using an 
invar material, for instance), and then measure the force at two accessible temperatures 
in the laboratory, for instance 300 K and 350 K. It would be even better if experiments 
could be carried out at liquid He temperatures; we understand such experiments are 
underway. 
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